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We study the spectral properties of thermal fluctuations on simple liquid surfaces, sometimes 
called ripplons. Analytical properties of the spectral function are investigated and are shown to 
be composed of regions with simple analytic behavior with respect to the frequency or the wave 
number. The derived expressions are compared to spectral measurements performed orders of 
magnitude below shot noise levels, which is achieved using a novel noise reduction method. The 
agreement between the theory of thermal surface fluctuations and the experiment is found to be 
excellent. 

I. INTRODUCTION 

Thermal fluctuations are ubiquitous and exist for practically everything we see and touch. However, they 
tend to be too small to be observed or measured directly, except under special circumstances. Phenomena in 
which thermal fluctuations can be examined directly in non-exotic materials are surface fluctuations of liquids, 
sometimes called "ripplons" [l| . Using surface light scattering, thermal surface fluctuations of liquids have been 
studied for some time[2|. Other direct thermal fluctuation measurements include high power interferometry of 
mirror surfaces Q and fluctuations of surfaces with exceptionally low surface tensionQ. 

In this work, we detect reflected light from surfaces to measure their inclinations^- 7] and apply this to 
thermal surface fluctuations of simple liquids. The measurements are somewhat complementary to the spectral 
measurements performed at specific wavelengths. The measurements require less power, can be performed on 
smaller samples and can be applied to strongly viscous fluids. Furthermore, the measurement system can be 
simpler. While shot noise is often regarded as an unavoidable limitation in such measurements, we show that 
such is not the case. The measurements involve novel methods that allow us to measure the spectrum directly, 
down to several orders of magnitude below the shot noise level. As we explain, this principle is not limited to 
thermal noise measurements nor to surface light scattering measurements. The experimental results agree with 
the theory quite well. In the process, we elucidate the simple analytic behavior of the spectra and show how it 
is reflected in the full spectrum, which can be seen in experiments. 

The dispersion relation for surface waves on a simple liquid surface can derived from the Navicr-Stokcs 
equation, when the viscosity of the liquid can be ignored, as 



uj(k) = ^^+gk . (1) 

Here, g is the gravitational acceleration, p, a are the density and the surface tension of the liquid, k is the wave 
number and uj is the (angular) frequency. The viscosity becomes more important for shorter wavelengths and 
dissipation will play an essential role below. Gravitational effects are more important for longer wavelengths 
and is negligible for wavelengths much smaller than it a / (pg) . For liquids we examine, namely, water, ethanol 
and oil, gravitational effects are unimportant for scales below 10 mm. The samples we examine have surface 
sizes of few mm and longer wavelengths are effectively cut off, so that henceforth we ignore effects due to gravity. 
Under these circumstances, the dispersion relation reduces to the following two equivalent relationships. 

pur x 




uj = u) R (k) = \ , k = k R (uj) = . (2) 



In what follows, we examine the full spectrum both experimentally and theoretically, including the full effects 
of dissipation. 

The paper is organized as follows: In Sect. UH we explain the surface light reflection experiment and derive 
what precisely is measured by this. The properties of thermal surface fluctuations of simple liquids are examined 
in Sect. HTT1 and the approximate simple analytic behavior of the spectra are derived. The shot noise level in our 
experiments are assessed in Sect. II VI The noise reduction method explained in Sect. fVl allows us to detect weak 
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signals buried under the shot noise level. The general applicability of the principle is also clarified. Finally, we 
combine the theoretical and the experimental results in Sect. IVII and find that they agree. The limitations in 
the experiment and further directions for research are also discussed. 



II. THE EXPERIMENT AND THE MEASUREMENT 



In the experiment (Fig.[TJ, inclination fluctuations of a liquid surface are measured using two sets of measure- 
ment systems each with a laser beam 1 (wavelength 638 nm) and 2 (658 nm). The average inclination within the 
beam area essentially acts as an optical lever and the inclination deflects the laser beams to be detected by the 
dual-element photodiodes (DEPD 1,2, S4204 Hamamatsu Photonics, Japan) ;5]. The two sets of photodiodes 
are necessary here to eliminate the random noise statistically, as explained in Sect. [Vj The DEPD signals are 
amplified and then fed into a computer via analog to digital converters (ADC, 14 bit, ADXII14-80M, Saya, 
Japan). Fourier transforms and averagings are performed by the computer. 

The laser beam power at the sample is 0.5 mW each. The beam is focused down to order and its diameter 
(~ X/NA) is varied by changing the numerical aperture (NA) of the objective lens, where A is the wavelength 
of the probe laser beam. This beam size is considerably smaller than that used in the standard light scattering 
methods. The amplitude of the waves are small compared to the wavelength so that the reflected light is almost 
all collected by the objective lens. Therefore, compared to the standard light scattering experiments which 
observes only a small fraction of the scattered light, we obtain larger signals for a given beam power. This can 
be crucial in low power measurements. 
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FIG. 1: Two laser beams with wavelengths 638 nm and 658 nm are combined at a dichroic mirror (DM1) and focused by 
a microscope objective lens onto the liquid surface. The reflected light with different wavelengths are separated at DM2 
to obtain two independent inclination measurements of the same surface. The inclination is converted into electric signals 
using DEPDl, DEPD2 and are fed via ADC's into a computer where FFT (fast Fourier transform) and averagings are 
performed. A polarizing beam splitter (PBS) and a quarter wave plate (QWP) are included to extract the light reflected 
back from the sample efficiently. DM3 is used for viewing the sample through a video camera. 



We now explain what is measured in this experiment and how this relates to the spectrum of surface fluctua- 
tions. We measure fluctuations in the average inclination of the surface. Denoting the surface level as <fi(r, t), the 
average inclination a±(t) can be obtained by approximating <f>(r,t) by a linear profile, minimizing the following 
quantity: 

r d 2 rG{r)\<p{r,t)~{a Q {t)+ ai {t)x)\ 2 . (3) 

Here, r — (x, y) are coordinates on the surface and G(r) is the beam intensity profile. This leads to the 
expression for the average displacement and inclination, 

^ -l 

2„ 



oo(t) = C / d 2 rG(r)<l>(r,t) , C =[ d 2 rG(r) , (4) 



oi (t) = Ci J d 2 rG(r)x(/)(r,t) , C x = (J d 2 rx 2 G(r)^J 



(5) 



While we can solve for a,j(t) more generally, we assumed here that the profile is symmetric, J d 2 rxG(r) = 0, 
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which applies to our case, as discussed below. We define the Fourier transform of a,j as 

a» = -J= f T/2 dte^a 3 {t) , j = 1,2 , (6) 

Vi J-T/2 

where T is the measurement time, which is much longer when compared to the other time scales involved. 

Since the correlation function of the surface fluctuations is translation invariant both in space and time, it 
can be expressed as 

<*(r, t)<t>(r', f )) =J^J-^2 e^-r)-M t -t') P{K u) . (7) 

Here, P(k,uj) is the spectral function of the fluctuations of the surface displacement and (• • • ) denotes the 
statistical average. Using this correlation function, the fluctuations in the inclinations are obtained as 

P(k,Lu) . (8) 

Since we are using a laser beam well described by a Gaussian profile to observe surface fluctuations, the profile 
function G(r) can be expressed using the beam diameter b as 

G(r) - I Q e- Sr2 / b2 , (9) 

The spectrum measured in the experiment is (we henceforth use the notation ui = 2nf) 

S(f)= 47r<|ax (a,)| 2 ) , (10) 

taking into account that the measurement is in frequency space and the one-sidedness of the spectrum. Combin- 
ing the results above, we derive a compact expression for the fluctuation spectrum observed in the experiment, 

POO 

S(f)= / dkk 3 e- b2k2 / 16 P(k,27rf) . (11) 
Jo 

It should be noted that this formula applies to general surface fluctuation spectra measured using this method 
and is not limited to liquids nor to thermal fluctuations. This result specifies the measured spectrum S(f) com- 
pletely including its magnitude, given the spectral function P(k, uj) and the beam diameter b, and is independent 
of the beam power applied. As can be seen from the expression, the role of the beam size is to effectively cut off 
the k integral of the spectral functions for values over ~ 27r/6. This occurs because the inclination is effectively 
averaged within the beam spot, so that shorter wavelengths are effectively averaged out. It also explains why 
we integrate up to infinity in this formula; while, in principle, the wavelengths of surface fluctuations should 
be cutoff at atomic length scales, this is much smaller than b so that using infinity as the upper limit in the 
integration region introduces negligible difference, due to the Gaussian damping. The lower limit of the inte- 
gration region should, strictly speaking, be set to ~ 2tt/L, where L is the size of the sample (L is few mm in 
our experiment) providing an upper bound for wavelengths. However, the difference from setting the lower end 
of the integral to zero as in the formula can also be ignored, as will become clear below. 
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III. ANALYTIC STRUCTURE OF THE SPECTRUM 



The spectral function for thermal fluctuations of simple liquid surfaces has been derived previously [. 
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While this expression for the spectrum is analytic, its behavior is not apparent. The behavior can be split into 
several regimes depending on the viscosity. We summarize the behavior concisely below to gain insight into the 
behavior of the spectral function and for later use. 



A. Leading analytic behavior with respect to k 



Let us study the behavior of the spectral function P(k, w) with respect to k with uj fixed. The dimensionlcss 
measure of viscosity of the liquid, rfui/(po~) influences the behavior of the spectral function P(k, uj) qualitatively. 
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Any liquid is highly dissipative for high enough frequencies, which is intuitively natural. When the viscosity is 
effectively low, rfuj/(pa) < l/(8-\/2), the leading order analytic behavior P (k,oj) can be obtained as 



k n T 
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(13) 



k > 2- 1 / 6 k R (u) 



In this case, P(k,ui) has a peak close to k = k R (oj), which becomes less prominent as T] 3 uj/(pa) increases. When 
r] 3 uj/(pa 2 ) > 1/(8-^/2), the leading behavior of the dispersion relation splits into three regions as 
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In both cases with low and high viscosity, for long wavelengths compared to 2ir/k R (uj), the spectral function is 
governed by viscous behavior and depends on p. On the other hand, for short wavelengths, it is suppressed by 
the surface tension of the liquid and is independent of its density. 

In Fig. [2] we compare the full spectral function Eq. (1121) with its approximate analytic behavior derived above. 
We see that in all the cases, the spectrum is well reproduced by the simple analytic behaviors, except for the 
peak seen in the water surface fluctuations at lower frequencies. The peak behavior reflects long lived waves 
and disappears when the viscosity is effectively high, due to dissipation. In both Eqs. (fT3"|) . (fT4")l . the spectral 
function is independent of / for large k, which can be seen from the plots. The simple analytic formulas derived 
above capture the situations with weak and strong viscosity, which have qualitatively different behaviors. 
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FIG. 2: (Color online) Behavior of P(k, 2nf) for (a) water and (b) oil with respect to k at fixed /. Frequencies are 
/ = 10 3 , 10 5 , 10 7 [s _1 ] (red, green, blue lines respectively), with the maximum of the spectrum being larger for higher /. 



The corresponding simple analytic behaviors Eqs. (I13p . (|14p are also shown (thin lines, magenta) and matches well with 
the full spectral function and are almost invisible, except at the boundaries between the regions. 



B. Leading analytic behavior with respect to uj 

We now analyze the behavior of P(k, to) with respect to lu for fixed k. For a liquid with low viscosity effectively, 
rj 2 k/(pa) < 1/(4 v^), the leading order analytic behavior Po(k,u)) is broken up into two regions as 
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For the highly viscous case, the leading analytic behavior can be broken down into three regions thus. 
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FIG. 3: (Color online) Behavior of P(k, 2nf) for (a) water and (b) oil with respect to frequency / at k = 10 3 , 10 s , 10 7 [m _1 ] 
(red, green, blue lines respectively). The spectra with larger k have more fluctuations at higher /. Their analytic 
approximations Eqs. (|15[1 . (|16p (thin lines, magenta) agree well with the full spectral function and are almost invisible, 
except at the boundaries between the regions. 



C. Analytic behavior of the integrated spectrum 



While the spectrum S(f) in Eq. (jllj) can be computed numerically, we now clarify its rough analytic behavior 
using Po(k,2irf) derived above. In the integrated spectrum S(f), b provides a cutoff fc max for the integral, 
whose relation we specify below. The spectrum can be approximated by 
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This can be computed explicitly. For the low viscosity case rfuj/(pa 2 ) < 1/(8v2)j 
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Taking into account the formula for S(f) in Eq. (fTTj) . we use fc max = 2 1 / 4 4/6, which satisfies 
J °° dk exp(— b 2 k 2 /16)k 3 = J Q fcmax dk k 3 . In Fig. 21 the spectrum S(f) is compared to its approximate ana- 
lytic behavior So(f) in Eqs. (|18p . fp~9|) . It can be seen that the essential features of S(f) are well reproduced by 
the simple analytic formulas. Examining in more detail, we see that the agreement is better for more viscous 
liquids. This is presumably due to the peak contribution unaccounted for in the formulas Eqs. (fT3|) . (fl4l) . which 
do not exist for more viscous fluids. It can be seen from Eqs. (|18[) . (fT9|) that S(f) behaves as ~ rj/a 2 for lower 
frequencies explaining why oil has larger fluctuations than water in this regime. Also, from the formulas, we 
can see why the fluctuations are larger for smaller b (or equivalently, larger fc max ) over the whole spectrum, 
both when the viscosity is weak and strong. 

We note that the calculations above also show why the integration region in the spectrum formula Eq. (jlip 
can be taken down to zero as long as the cutoff is well below the upper limit of the region, since had we put in 
a lower cutoff fci ow , its contribution to S(f) would behave as ~ kJ OVJ . This condition is equivalent to the surface 
size being much larger than the beam diameter, which is always satisfied in our experiments. 
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FIG. 4: (Color online) Behavior of the integrated spectra of surface fluctuations S(f) for (a) water and (b) oil with 
respect to /. The spectra are shown for 6 = 1, 2.5, 10 [/im^ 1 ] (red, green, blue lines respectively), along with their 
analytic approximations (thin lines, magenta). The spectra for smaller b have larger fluctuations across the spectrum. 



IV. SHOT NOISE LEVEL 



While our measurement precision is not limited by the shot noise of the probe laser beam, it is of import to 
understand how the shot noise level would appear in our setting. The electric current through the photodiode 
is / = neP/his, where V is the signal power, e is the electron charge magnitude and k is the quantum efficiency 
of the photodiode. k ~ 0.8 for the photodiodes we use. The current generated by shot noise is 



7 SN = v/^TA/ ^e^j^-Af . (20) 

A/ is the frequency range of the measurement. The signal in the experiment comes from difference in the 
current AI due to the geometric effects of the optical lever. For an inclination 9, the signal is 

AI = K - . (21) 

hvNA v ' 



The shot noise level in our experiment is the size of the angular fluctuations 0| N corresponding to the shot noise 
current. This is what appears in the measurements, had we not used the noise reduction through correlations 



described in Sect. [Vj 6*g N can be obtained from AI ~ Tsn to be 



0§n - -^r A f ■ (22) 

Here, the current I is the photoelectric current collected by the dual photodiodes. 6>g N is smaller for larger 
signals, which is natural. Also, for larger numerical apertures, the beam spot is focused down to a smaller 
size so that the effect of the optical lever is smaller, hence the shot noise effects are larger. In our experiment, 
I ~ 0.5 fj,A so that 6>g N ~ NA 2 x 1.6 x 10~ 13 Hz^ 1 . It is crucial to separate out the signal from this noise, using 
methods explained in the next section. 
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V. NOISE REDUCTION THROUGH CORRELATIONS 

An essential feature of thermal fluctuations is that they are random. Since our objective is to directly measure 
the fluctuations under normal circumstances, the fluctuations are furthermore small. Even in ideally executed 
experiments, some random noise, such as shot noise, always exist. Therefore, to measure weak random signals, 
we need to separate out the random signal from the random noise. This is possible under rather general 
circumstances, as we now explain [7J. The principle behind this noise reduction is not limited to thermal 
fluctuations or optical measurements, but applies generally to the extraction of random signals from random 
noise. The conditions for its applicability will be discussed below. 

A detector measurement Di = S + Ni consists of the desired signal S and some noise Ni, independent of S. 
Denoting Fourier transforms with tildes, the power spectrum obtained under simple averaging is 

(\D 1 \ 2 ) = {\S\ 2 ) + {\N 1 \ 2 ) . (23) 

Since the signal itself is random in nature, there is no way to distinguish the signal from the noise, so that the 
signal can not be measured, unless the signal is larger than the noise, (\S\ 2 ) ^> (|Aq| 2 ). If we use only one 
measurement, this is an essential limitation. 

To overcome this obstacle, we make another independent measurement of the same signal, D2 = S + N2 , 
where N2 denotes the noise for this measurement. Then, 

(IhDi) ^ (\S\ 2 ) , AA^co , (24) 

eliminating the random noise Ni,N%. Here Af is the number of averagings. This results holds since the 
measurements D\ , D2 are independent and the cross-terms of decorrelated random observables (and their Fourier 
transforms) vanish under averaging. The relative error in this method is ~ I/x/AT, which arises from the 
statistical nature of the method. In principle, given enough averagings, we can suppress shot noise and other 
random noise effects to an arbitrarily small size. 

The crucial requirements for our method to work is that multiple independent measurements can be made and 
that the signal is stable enough to withstand averagings. The independence of the measurements, or equivalently, 
the decorrelation of the noise in them is clearly crucial. As the signal becomes weaker, stricter independence is 
required, which in practice can be quite delicate. For instance, cross talks can arise in electronic circuits, unless 
they are completely separated and electronic signals can affect each other through electromagnetic fields in the 
intervening space. These properties put practical limitations on the reduced noise level. While our method can 
not be used for an one-time event, it can be used for any recurrent signal. 

Another possible approach to reducing the relative noise is to increase the signal strength. In our context, 
this would mean increasing the beam intensity. However, this is not always applicable, since a stronger beam 
will affect the sample. Even for simple liquid surfaces, it leads to more evaporation and can lead to less 
precise measurements. More generally, if we consider biological materials or medical applications [9], using a 
strong light source is often excluded. Correlation measurements have been used previously in surface light 
scattering experiments 10]. Our approach differs from those in that we use the cross-correlation of independent 
measurements of the same signal to reduce the noise. 

VI. EXPERIMENTAL RESULTS AND THEORY 

In Fig. [5l we compare the experimental results with the theory explained above for water, ethanol 
and oil with various beam sizes. The fluid properties we used for water, ethanol and oil are 
(p[kg/m 3 ],cr[kg/s 2 ],r?[kg/(m-s)]) = (1.0 x 10 3 ,7.3 x 10~ 2 , 1.0 x 10~ 3 ), (0.79 x 10 3 ,2.2 x 10~ 2 , 1.1 x 10~ 3 ) 
and (0.92 x 10 3 ,3.0 x 10~ 2 , 0.124), respectively. The agreement between the theoretical formula Eq. (1 1 1 [) and 
the experimental measurements is mostly quite satisfactory, including its beam size dependence. There is some 
excess signal in the water surface fluctuations at low frequencies for small b (b=1.3 /im case in Fig.(5ja)), whose 
cause is explained below. The rough features of the spectra can be understood from Eqs. (fT5|) and (|T9l : the 
spectrum behaves as S(f) ~ i]/(a 2 b) at low frequencies. So, the fluctuations are the largest for oil due to its 
large viscosity and smallest for water due to its large surface tension. The fluctuations are larger for smaller b. 
At high frequencies, as we decrease b, the effective cutoff for the wavelength decreases and fluctuations at larger 
frequencies become more apparent. 

The excess measurements in the water surface fluctuation spectrum at low frequencies for b — 1.3 /jm can be 
explained as follows: Ideally, the experiment measures only inclination fluctuations of the surface. However, 
given the high sensitivity of the measurements, vertical displacements can mimic inclinations due to the focusing 
of the objective lens, when the beam alignment is not perfect. This effect is contained in (|a (w)| 2 ) which behaves 
as ~ 1/f at low frequencies. For smaller b, this effect is larger due to the steeper focusing and can be seen in 
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the water surface fluctuation measurements at low frequencies, since for a given beam size, S(f) is smallest for 
water. 

The shot noise levels in the measurements are 6 x 10~ 14 , 3 x 1CP 14 , 2 x 10~ 15 [Hz -1 ] for the beam diameters 
b = 1.5, 2.5, 10 |/im], as explained in Sect.|lV] Therefore, we see that the noise reduction using signal correlations 
explained in the previous section is crucial for examining even the qualitative features of the spectra, since most, 
and in some cases all, of the spectra is below shot noise levels. 

Some comments regarding the calibration of the measurements is in order. On the theoretical side, given 
the properties of the liquid and the beam diameter, there are no further parameters at all in the spectrum 
Eq. (fTTj) and is specified completely. Experimentally, the frequency dependence of the spectrum can be measured 
precisely. Calibrating the overall magnitude of the measured spectrum is more difficult and this is done using 
a piezoelectrically driven mirror with a known oscillation amplitude. While this works well when b is large, for 
smaller b, the shallowness of the depth of field makes the calibration less accurate. Another complication is that 
the liquid in general evaporates during the measurement so that the beam defocuses. This, in effect, increases 
the beam diameter and can influence the spectral shape. This problem is clearly more acute for larger beam 
powers. A typical measurement of simple liquid surface fluctuations takes around 20 seconds and the beam 
power applied is 0.5 mW. Given the excellent agreement of the observed fluctuation spectra of simple liquid 
surfaces and theory, it might be reasonable to use thermal surface fluctuations of a well studied specific liquid 
such as oil to finely calibrate the system when applying the measurements to more general samples. 
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FIG. 5: (Color online) Observed fluctuation spectra for (a) water (b = 1.3, 3.1, 11.5 [/an]), (b) ethanol (b = 
0.9, 2.5, 9.5 [/im]) and (c) oil (b = 1.7, 2.5, 11.0 [Atm]), (red, green, blue lines from smaller to larger 6). Respective theoret- 
ical spectra S(f) are also shown (thin lines, magenta). The spectra for smaller b have larger fluctuations. 

In this work, we directly studied the spectra of thermal fluctuations for simple liquids, using surface light 
reflection methods. The spectra obtained are integrated over wavelengths and we have also investigated the 
dependence of the spectra on the beam diameter. In the process, we applied a novel general method for noise 
reduction, using the correlation of independent measurements of the same signal. The spectra obtained matches 
well with theory, whose analytic behavior can be summarized rather simply. The measurement method is 
complementary to the surface fluctuation measurements for specific wavelengths, when applied to simple liquids. 
The method, especially when combined with the noise reduction method, has a broad range of applicability 



9 



since the required sample size, observation time and power are small. 



[1] M. von Schmoluchowski, Ann Physik 25, 225 (1908); L. Mandelstam, Ann. Physik 41, 609 (1913). 

[2] D. Langevin, "Light scattering by liquid surfaces and complementary techniques", Marcel Dekker, New York (1992). 

[3] K. Numata, M. Ando, K. Yamamoto, S. Otsuka, K. Tsubono, Phys. Rev. Lett. 91, 260602 (2003); E.D. Black et al, 

Phys. Lett. A 328, 1 (2004) 
[4] H.J. Lauter et al, Phys. Rev. Lett 68, 2484 (1992); A. Vailati, M. Giglio, Nature 390, 262 (1997) 
[5] T. Mitsui, Jpn. J. Appl. Phys., 47, 6563 (2008) 
[6] A.Tay et al, Rev. Sci. Instrum. 79, 103107 (2008) 
[7] T. Mitsui, K. Aoki, Phys. Rev. D E 80, 020602(R)-1— 4 (2009). 
[8] M.-A. Bouchiat, and J. Meunier, J. de Phys. 32, 561 (1971). 
[9] W. Denk and W. W. Webb, Appl. Opt. 29, 2382 (1990) 
[10] For instance, J.C. Earnshaw, Appl. Opt. 36, 7583 (1997). 



